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Abstract 

We investigate the stabiHty of the extremal black p-brane which contains a n-form and 
a dilaton. We show that the instability due to the s-mode, which was present in the 
uncharged and non-extremal p-brane, disappears in the extreme case. This is shown to be 
consistent with an entropy argument which shows that the zero entropy of the extremal 
black hole is approached more rapidly than the zero entropy of the black p-brane, which 
would mean an instability would violate the second law of thermodynamics. 
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1. Introductory remarks. 

In four dimensions, black holes are characterised by a small set of parameters: charge, 
mass and angular momentum. The simplest solution, Schwarzschild, representing an un- 
charged non-rotating black hole has long been known to be classically stable, and in Ein- 
stein gravity, the stability of the exterior spacetime of the charged and rotating coun- 
terparts (Reissner-Nordstrom, Kerr-Newman) has also been demonstrated. In this latter 
case, there is an instability of the inner Cauchy horizon, however, this is hidden from us 
by the event horizon, and is therefore presumably not of physical interest to the exterior 
observer. Therefore, once a black hole is formed, it remains, forever shielding the details 
of its constituents from the outside observer, who may only measure the charge, mass, and 
angular momentum of that matter which formed the black hole. 

Quantum mechanically, we have a very different picture. Black holes become black 
bodies, semi-classically radiating a purely thermal spectrum. In theory, a black hole could 
radiate away completely, possibly losing forever all the detailed information of its con- 
stituents. Such a picture, in which quantum mechanical unitarity is lost, is potentially 
unsettling, and certainly begs the question as to whether there is not some modification 
of semi-classical gravity in which information retrieval is in some way possible. An ex- 
amination of black holes in low energy string theory is a step in this direction. Here, 
the Einstein action is modified by the addition of a dilaton with a non-trivial coupling to 
'electromagnetism' 



Such 'dilatonic' black holes have several attractive features. In Einstein gravity, charged 
black holes have an instability of the inner Cauchy horizon due to matter perturbations in 
the exterior spacetime^, however, there is no static charged black hole solution in Einstein 
gravity with only one horizon and a spacelike singularity. On the other hand, in low 




(1.1) 
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energy string gravity the dilaton greatly changes the causal structure of charged black holes 
making them like Schwarzschild, having one event horizon and a spacelike singularity^. 
This structure is generic, even if the dilaton has a mass^, as is generally believed from, 
although not necessarily implied by, the principle of equivalence. Additionally, extremal 
magnetically charged dilaton black holes have an infinite throat, which, unlike the spatially 
infinite throat of extremal Reissner-Nordstrom, is a true internal infinity as seen in the 
so-called string metric 

^string ^ ^Einstein (l'^) 

Such an internal infinity allows an eff'ective dimensional reduction down to two-dimensions 
in the throat and has been a topic of much debate as a simplified model for black hole 
accretion / radiation. 

Since superstring theory is most naturally, though not necessarily, associated with a 
ten-dimensional background, ideally one should examine string gravity in higher dimen- 
sions. In four dimensions, an event horizon must be topologically spherical, but in higher 
dimensions this need not be the case, for example, in ten dimensions we can have S'^x 
H^, S^x WiP etc., topologies for the horizon. In previous work^'^, we have pointed out 
that such black p-branes are unstable for a large range of (magnetic) chargermass ratios, 
however, we were unable to demonstrate an instability for the extremal black p-branes. 
This was largely because our numerical analysis was tailored towards a general charge:mass 
ratio, and if "Q=M", a degeneracy occurs in the perturbation equations reducing the di- 
mensionality of the numerical problem from 5 to 3, this degeneracy causing numerical 
instability in the solution of the equations. In this paper, we seek to close this loophole 
in our argument, tailoring our analysis specifically to extremal black p-branes, providing 
a convincing argument for the stability of such objects. 

The layout of the paper is as follows. We first review some heuristic arguments, 

3 



particularly a thermodynamic entropy argument, as to the stability (or otherwise) of black 
p-branes. We then set up the formalism involved in solving the problem and argue that 
the only possible instability is an s-mode. We then examine this s-mode numerically and 
show there is no solution to the perturbation equations corresponding to an unstable mode. 
Finally we summarise our results and discuss their implications. 

We would first like to examine several arguments in favour of stability for extremal 
black p-branes. It is generally believed that, being supersymmetric and satisfying an appro- 
priate Bogomolnyi bound, extremal black p-branes should be stable, although formulating 
an appropriate argument in low energy string gravity would be necessary to prove this. 
One of the main problems is that the mass appearing in the Bogomolnyi bound is generally 
the mass per unit p-area, and this is defined by taking the ADM mass of the D-dimensional 
black hole. Such a procedure relies crucially on the translational invariance of the metric 
in the p orthogonal directions. Since our instability, when it exists, involves a dependence 
on these directions, it technically sidesteps the bound. In order to close this route one 
would have to define a generalised notion of 'asymptotically fiat', re-derive the ADM mass 
and positivity theorems, and finally the Bogomolnyi bound for 'approximately p-planar' 
objects. 

A second consideration is that certain extremal p-branes are relevant in the context 
of S-duality - a generalised strong-weak coupling duality which, loosely speaking, swaps 
solitons and gauge particles, electric and magnetic charges. If a supersymmetric brane 
were unstable, it would be rather embarrassing for the dual particle! Although this is 
not a good reason for stability, there is sufficient circumstantial evidence in favour of a 
limited form of S-duality that it would be most remarkable if we were to discover that the 
supersymmetric branes were unstable. 

Finally, our main reason for expecting stability is a thermodynamic one. One of the 
most convincing arguments that we should expect an instability of the uncharged black 
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p-brane was an entropic inequality. We calculated the entropy of a length L of black 
string in five dimensions and compared it to that of a hyperspherical black hole. For the 
string, the entropy was proportional to M^L, for the black hole, M^/^L^/^. Regarding M, 
the mass per unit length of the string as a constant, this showed that the entropy of the 
hole is greater for L > 0{M), and hence that a series of black holes is preferred to the 
string. However, such an argument did not include charge and the non-trivial dilaton field 
that this entails. It is therefore worth re-investigating this argument and generalising it 
to include charge. We will perform this analysis for general N, the number of spacetime 
dimensions, and general D = N—p, so asto be able to show the genericity of the argument. 

In order to calculate the entropy of the black p-brane, it is most convenient to work 
in the "Einstein" frame, i.e., to make a conformal transformation so that the gravitational 
part of the action is in Einstein form 



d'^ x^/^R + .... (1.3) 



The appropriate conformal factor for N dimensions is 

^ -40/(iV-2) 

t7Einstein ^/string V * / 

where cj) is the dilaton field. Applying this conformal factor to the metric of a p-brane 
yields: 



dst = ^ ^dt' + 



Bin 



(1 _ (!-)I^-3)l^ (1 - (Ii)^-3)(1 - (^)^-3)^ 



(1.5) 



+ r2(l - {•—)^~^)^- d^l_^ + (1 - i^)''-^)^- dxM 
^From this, we read off the entropy using S — \A as 

S=\An-.r^^--{l-{^r-')^' (1.6) 



where 

Tl+l 

An = (1.7) 

r(^) ^ ^ 

is the area of a unit n-sphere. 

In order to compare the entropies of a p=(N-D)-brane with an N-dimensional black hole, 
we need the entropy as a function of mass and charge. To find the mass of a D-dimensional 
black hole, note that as r — > oo 



IGtt Md 
{D - 2)Ad-2 H 



l^ool ~ 1 - TT^^TVl (1-8) 



Reading off the asymptotic form from (1.5) gives 

(D - 2)Ad-2 



M 



D 



IGtt 

Additionally, Q'jj = ^^^(r+r_)^~^, which gives 



D-3 _ Ar-4 -D-3 
^+ N-2'- 



(1.9) 



_ SttMd / 6A7r^Ml 2Ql{N - A) 

(D - 2)Ao-2 V - 2)^^1.-2 iN-2){D- 3) 



(1.10) 



D-3_ -8ttMd{N- 2) / 647r2M2(7V-2)2 2Ql{N - 2) 



(D-2)(Ar-4)ylz,_2 V (^-2)'(^^-4)2^1)-2 (A^-4)(D-3) 



Thus we input the appropriate form of r± into S'd to find the entropy per unit area of 
the black p-brane, and setting D = N gives the relevant quantities for the N-dimensional 
black hole. We then compare the entropies by taking a 'cube' of p-brane, side-length L, 
and setting 

Mn = L^~^Md 

(1.11) 



n - / (^-3) (i^-2) An_ON-D^ _ /^TTT^ 
- V(i^-3) {N-2)A^_2^ - ^2^^ - 



D 



iJV-2 



The rather unusual factors appearing in the charge formula are to ensure that the extremal 
limit of the black p-brane corresponds to the extremal limit of the black hole. If this were 
not the case, then the p-brane would either become stable before reaching its extremal 
limit (a scenario we have already eliminated), or would be unstable at all length scales at 
its extremal limit, and in particular at small length scales close to the extremal limit; this 
too is ruled out by our previous work which showed that the length scale of the instability 
was tending to infinity as the extremal limit was approached. We therefore assume that 
these limits coincide, and introduce the constant c as representing a general charge:mass 
ratio, so that c — > 1 as the extremal limit is approached. 

We now set L^~^Sd = to find the critical length, L*, at which it becomes 
entropically favoured to form an N-dimensional black hole. Taking 2tt/L^ gives the 
critical frequency which, after some algebra, is found to be: 

1 -(Ar-3) 

(D-2)(iV-3) / r-l N-l\ \ N-D / r( D-l\ \ {D-3)(N-D) , , , ^ 

,, _ 277 (N-^\ In-2Mn-^ [ ^\ — ) \ i ^\ — ) \ ( 2{D-2) \ N-D 

^*~^^^ 2 j y N-2 J y D-2 J {(N-2){N-4)) 

-(ND-2N-3D + 7) jV-3 

X (l+^l+c2(Ar-2)(Ar-4)) {N-2)(D-3) (jV-3-^l+c2(iV-2)(Ar-4)) 

(1.12) 

Clearly 

/X* oc(A^ - 3 - Vl + c2(iV-2)(iV-4))^ 

= (AT - 3 - Vl - c2 + c2(Ar - 3)2)^ (1-13) 

— > as c ^ 1. 

Figure 1. shows the numerical and analytical value for Umax for which the instability just 
disappears. It is seen that in both cases Umax goes to zero in the limit of extreme charge. 
Hence, as the extremal limit is approached, the scale at which it becomes entropically 
favourable to form a black hole tends to infinity. Thus the entropy argument, which was 
so convincing for uncharged black holes, actually reinforces the idea that extremal black 
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holes should be stable. We now turn to setting up the formalism for the stability analysis. 
Since this has been done in detail in [5] , we will merely paraphrase the process here, quoting 
only the main formulae. 

2. The stability analysis. 

The low energy string action used by Horowitz and Strominger^ is 

J d''xV^e-'^[R + 4(V0)2 - ^^^F'] (2.1) 

where F is a (-D — 2)-form field strength, and (p the dilaton. This has an extremal mag- 
netically charged black p-brane solution: 

e-^-^ = (1 - (^)^-3) (2.2) 
r 



F=\l ^r^-'eo-2 



where the index i runs from L> + 1 to 10 and e£)_2 is the area form of a unit (D-2)-sphere. 
In order to write the perturbation equations, we use the usual notation 

Sgab = Kb , (2.3) 
and a conventional gauge choice of transversality 

^ah''^ = ^aih^^ - ihg'^'') = . (2.4) 

To perform the stability analysis, we Fourier decompose the perturbations in terms of 
the symmetries of the background spacetime. We perform our analysis in Schwarzschild 
coordinates, and transform to the tortoise coordinate at the future event horizon to check 



the regularity of the perturbation equation (see [5] for a more detailed discussion of the issue 
of boundary conditions and the generalised tortoise coordinate). The Fourier modes in the 
time and p-brane coordinates are of the form 6^*+*^*^^' for an instability. The spherical 
harmonic modes will depend on the number of dimensions, D, that the black hole sits in, 
as well as the tensorial nature of the perturbation we are analysing. However, we note that 
in the uncharged black p-brane, the higher angular momentum mode equations have no 
unstable solutions, and also that the s-wave instability shrinks in parameter range as we 
increase the charge (see figure 1.). In other words, adding charge has a stabilising influence. 
Thus we do not expect higher angular momentum perturbations of charged black holes to 
exhibit instabilities. Therefore, for the extremal black p-brane, we focus on the s-mode 
which caused an instability in the lesser-charged branes. This mode takes the form 



= e 



fir) 



SF = 



H'\r) H^\r) H'''{r) 

H^^lr) H^\r) H^^'lr) 

H^j(r) iy^*(r) i^^^Jr) 











K{r) 











0.. 



K(r)/sin^6' 0. 



(2.5a) 
(2.56) 

(2.5c) 



where we have used the non-trivial result from [5] that 6F — 0. 

Using the gauge conditions, and (2.5b), the perturbation equations in tensorial nota- 
tion are 

D5(j) - h^^VaVbCp + 2/l"Va(/)V6(/) - 4Va5</)V"0 



I a2---(in—2 



^jj _ 4)1 ca2---aD-2^ d 

cd 







4 



_},cd rp 

l" acaz...an-2^ bd 



03...0u_2 



= 



(2.6a) 



(2.66) 
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Inputting the mode decomposition (2.5) yields the full set of equations which were 
written down in appendix B of [5]. We therefore set r_|_ = r_ in these equations, and 
without loss of generality we set D=4, since we demonstrated that the existence of an 
instability was not dimensionally dependent. In this case, the equations of [5] simplify 
neatly. An important feature is that the equations for /, /i^^ become independent of all 
the other variables. It is in this sense that the extremal limit is degenerate, the number 
of independent physical variables in the perturbation problem reducing from five to three. 
We can therefore look first at these functions to see if the instability persists. We should 
also point out that in the extreme case the behaviour of the modes near r_|_ is diflFerent 
than in the non-extreme case, this is because of the effect of the equality of r_ and We 
do the analysis for Z) = 4 in 5 dimensions but it can be generalised for higher dimensions 
and D. The equations for these modes are: 



(r - r+)2 (2r + r+)(r-r+) , (r - 2r+)r+ 

-/ / -0. ^ 



iy^3 1^2 



r+(3r+ -2r),^^ 2 



2r^(r — T-i-)^ 



l^rr ^ ^ 0(2.7a) 



(r - r+f ^„ ^ (-3r + 2r+)(r-r+) ^^, ^ 2{r - r+)^ 



= 0(2.76) 



^ly^ iy^4 iy^3 iy^5 



^4 V 2 
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Using (2.7b,c) we can obtain 



ir-r+)^ ^rr>> _ (r-3r+)(r-r+) ^^^, _ 4{r - r+f ^, 

(2.8) 



where = Jl^ + We can see that these equations depend only on // through A. If we 
would find an unstable mode for a given value of A this would imply one for = and 
Vt equal to A. This would correspond to an instability for a 4-dimensional dilatonic black 
hole! However as will be shown below, no such instability was found. 

For general A, we cannot directly determine if there exists a regular unstable solution. 
We must therefore resort to numerical techniques. These involve finding the regular solu- 
tion space at infinity, which forms a three dimensional subspace of the total solution space. 
We then integrate in, adjusting the initial configuration to see if one matches to a regular 
solution near the horizon. In practise, as described in [5], we will search for evidence of 
such a match, rather than the match itself. 

As r — > oo we can expand the regular solution as 

/ ^ /ie^^ 

K Kie^"" (2.9) 

where fi.,Ki, H^'^ are arbitrary constants. In fact one of them can be set to 1 due to the 
linearity of the equation. 

In the limit as r — > r+ the behaviour of the field is rather different than in the non- 
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extreme case. We assume the form 



/ = /(r - T+r 

K = K{r-r+f (2.10) 

and find that 



/ =A+{r - r+)"^ +A_{r- r+)"i + B+{r - r+)"^ 
B_{r - r+)"2" + C+(r - r+)"^ + C-(r - r+)"3" 

2 2 

{at - iy+ («! - iy+ n 

r;^ 



^+ ^+ ^+ 



for which 



(2.11) 



± ^ 13/ 4A2r2 

af = a = /3^^-2 = --±-^l + (2.12a) 



a± = a = /3-l = 7- 2- ±Ar+ (2.126) 



1 , 1 



af = a = p -1 = ^-2 = --±-^l + 4A^rl (2.12c) 

where A^, A_, B^, C_|-, C_ are constants determined by the integrated function with 
large r hmit specified in equation (2.9). Once a solution has been found by integration we 
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can find A_, B+, B_,C+,C- from the functions f,K,h'^'^ and their derivatives. By- 
varying the initial conditions we vary the ratio A_/A^ etc, until these ratios go to zero. 
That would correspond to a regular solution. We have used a Newton-Raphson method 
to investigate this problem and we were not able to find any regular solution. Therefore 
we should take /, K, h'^'^ to zero. 

The next step is to look at the possibility of an instability due to the other modes. 
Taking into account the previous result the equation for h^^ becomes 

- "' fe"" * fe"" + + " '7 ^ = (2-13) 

For Ar_|_ > 1/4 it is possible to show analytically that there is no regular solution to this 
equation but for smaller values we have to resort to a numerical analysis similar to the one 
above. We did not find any regular solution to equation (2.13). The equation for the mode 
h^^ is similar to that for h*'^. We must therefore take both /i^^ and /i*^ to be zero. With 
this result we can finally show that and h^^ must also be zero as their equations 

are of the form 



^^h''" -g{r)h''' + A''h''r = (2.14) 



where g{r) is an unimportant function of r which can be computed from eq (B.lle) of [5]. 
These equations do not possess regular solutions. Thus we have shown that the irregular 
modes which were found in [4] and [5] do not have a counterpart for the extreme case. 
3. Conclusions. 

In this paper we have shown that the s-mode leading to the instability of the uncharged 
and charged p-branes does not extend to the extreme case. We have made a detailed 
search for such a mode and we were not able to find one. This absence of an instability 
confirms our expectations based on an extrapolation of the charged p-branes studied in 
[5] and on the entropy argument presented earlier. Black p-branes therefore behave in a 
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thermodynamically consistent manner. We should stress that we have not proved stabihty, 
for we have neither studied higher angular momentum modes, nor exhibited completeness 
of our basis. Nonetheless, the fact that the s-mode which gave an unexpected instability 
for non- and lesser-charged branes does not give an instability for the extreme case, and 
also that adding charge and considering higher angular momentum modes has a stabilising 
effect, we regard as compelling evidence for the overall stability of extremal black p-branes. 

Using this result, it is possible to argue that charge might prevent p-branes from 
breaking. When the instability sets in in the non-extreme case, the apparent horizon 
expands in some regions and shrinks in others, but the charge per unit length is unchanged, 
at least to linear order. Thus it seems possible that these regions will evolve towards 
extremality and would therefore be stabilised. In this case the endpoint would appear 
to be a rippled brane, and we are investigating the possibility of such static solutions. 
However, the chargeless case does not appear to have such a mechanism and we still 
believe that these p-branes will break and reveal, at least temporarily, their singularity. 
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Figure 1. The critical frequency at which the hyperspherical black hole is entropically preferred to the 
black string. The dotted line corresponds to the numerical analysis done in [5] and the full one to the analytical 
argument presented here. The apparent crossing of the two curves in the vicinity of T — = T-\- ia probably due 
to numerical sensitivity. In both cases the instability does not seem to extend to the extreme case (r_ = r_|_). 
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